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Abstract 

I We derive an expression for the non-Gaussian cosmic-microwave-background 

Q\ I (CMB) statistic if defined recently by Ferreira, Magueijo, and Gorski in terms 

of the slow-roll- inflation parameters e and r]. This result shows that a nonzero 
value of if in COBE would rule out single-field slow-roll inflation. A sharp 
change in the slope of the inflaton potential could increase the predicted value 
0> ; of If, but not significantly. This further suggests that it will be difficult to 

account for such a detection in multiple-field models in which density per- 
turbations are produced by quantum fluctuations in the scalar field driving 
inflation. An Appendix shows how to evaluate an integral that is needed in 
. our calculation as well as in more general calculations of CMB bispectra. 

I> 

^ ■ I. INTRODUCTION 

• Ferreira, Magueijo, and Gorski (FMG) have recently found evidence for a non- 

Gaussian distribution of cosmic-microwave-background (CMB) temperature fluctuations in 
2 ! the COBE data (as have several other groups The common lore is that this result 

I is inconsistent with the nearly Gaussian distribution of temperature fluctuations expected 
^ ' from inflation. Although the result may be due to foregrounds or some curious systematic 
effect 0-0, it is still worthwhile to state more precisely the implications for inflationary 
^ . models if the non-Gaussianity is indeed in the CMB. That is the purpose of this paper.Q 

Inflation predicts the distribution of primordial perturbations to be very nearly Gaussian, 
but self-interactions of the inflaton field should in fact produce at least tiny deviations from 
Gaussianity P-[T3[]. We derive here an analytic expression for FMG's non-Gaussian statistic 



/; given in terms of the usual slow-roll parameters e and rj in single-field slow-roll inflation 
models. Our results verify the common expectation that the detected value of If is too large 
(by at least five orders of magnitude!) to be consistent with slow-roll inflation. Motivated by 
evidence for a break in the power spectrum of the galaxy distribution [p!^ -[l6|, it is natural 
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to consider inflation models in which the slope of the inflaton potential has a discontinuity 



17| , p!8[| . We show that the increase of the predicted non-Gaussian signal is insigniflcant. We 
infer from this that such a large value of If should also be hard to come by in multiple-fleld 
models in which density perturbations are produced by quantum fluctuations in the fleld 



driving inflation. (Counter-examples include the models discussed in Refs. |l9l-[2T[].) 

We begin by reviewing in Section II the calculation of the large-angle CMB power spec- 
trum. We then move on in Section III to the calculation of the CMB bispectrum and present 
our result for If in Section IV. Section V considers an inflation model with a discontinuity in 
the slope of the inflaton potential. Section VI provides a discussion. An Appendix provides a 
recursive technique for evaluating an integral involving the product of three spherical Bessel 
functions that is needed for our calculation and will be needed for more general calculations 
of CMB bispectra. 



II. TWO-POINT CORRELATION FUNCTION AND POWER SPECTRUM 

We start by reviewing the calculation of the angular two-point correlation function and 
power spectrum of CMB temperature fluctuations. The photon temperature perturbation 
at a spacetime point (x, r) can be Fourier expanded: 



AT, 



,n) = j rf3fce^'^-"A(k,n,r), (1) 



where n is the direction of photon momentum. We always set the observation point to be 
at the origin (x = 0) and at the present epoch (r = tq), so we will not explicitly write these 
two variables in the following derivations. By Legendre expansion, we have 



AT 



rf3A;A(k,n)= /rf3fc^(-i)'(2/ + l)V;(k)A/(A;)P,(k-n), (2) 

1=0 



where ipi}^) is the initial gravitational-potential perturbation and Ai{k) is the photon transfer 
function. We have used the fact that the photon evolution equation is independent of the 
wavevector direction k. For a stationary random process, we have 

(^(ki)^(k2)) = Pi'\k)SD{ki + k2), (3) 

(2) 

where the amplitude (k) is the primordial power spectrum. For a scale-free primordial 
power spectrum, we have 

pf(fc)(xP-^ (4) 

where index n = 1 corresponds to a flat scale-invariant spectrum, which is close to those 
favored by generic inflationary models. The CMB temperature pattern may be written in a 
spherical-harmonic expansion with coefficients, 

f AT 
aim = J dhYi„i{n) —{h). (5) 
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The angular two-point correlation function is 

/AT AT \ 2/ + 1 

e(fii, fi2) = ( — (fii) — (n^)) = E C, PKfii ■ fi2). (6) 

where 

Ci = {47cf J edkPf{k)\Mk)\' (7) 

is the CMB power spectrum. We have used Eqs. @ and @) to obtain Eq. @, and the 
spherical-harmonic addition theorem, orthonormality of spherical harmonics, and 

{o.hmxO*l^m^ = ^hl2^rmm2 Cl, (8) 

to obtain Eq. (|^. 

III. THREE-POINT CORRELATION FUNCTION AND BISPECTRUM 

If the primordial random density fluctuation is non-Gaussian, then the three-point cor- 
relation function is in general non-vanishing. We then have 

(V^(ki)^(k2)^(k3)) = Pfiki, k2, ks)6D{k, + k2 + ks), (9) 

where (^1,^2,^3) is the spatial bispectrum of the gravitational potential. The angular 
three-point correlation function for the CMB can be written as 

/AT, ,AT, ,AT, A 
e(ni,n2,n3) = ( — (fii) — (112)^(113)) (10) 

imi 

(ni)y/2m2 (112)1^3^3 (113), 

k,mi 

where 

{aiimiai,r^2ai,m,) = (47r)3(-z)'l+^2+'3 J ^3^^ ^3^^ ^3^^ 

xr,*„^(ko>^.:™.(k2)n:™3(k3)5D(ki + k^ + k,) 

3(3 



X Pf {k,, k2, k,)Ai^ {k,)Ai^ (A:2) A,3 (^3) . 



This last equation is obtained by using Eqs. (H), (H), and the spherical-harmonic addition 
theorem. By using 

1 

Sn{k^ + k2 + ks) = ^ e*('^^+i'^+'^«)-rf3^, (12) 
e^k.x ^ 4nY.^'Ukx)Y.yiUm*M), (13) 
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and the Gaunt integral, 



fdnv Y Y (2^1 + 1) (2^2 + 1) (2/3 + 1) //i kh] k k h\ .... 



where the (...) is the Wigner 3j symbol, we obtain 

^1 k I3 



{O'hrmO'hmiO'hma) ~ [ jjii m2 ] "^'I'^/s' (1^) 



where 



^3 (2/1 + 1) (2/2 + 1) (2/3 + 1) k k h 



B,,, - iS.r f - - - - [ - - - j (16) 

X kf dki kj dk2 kl dks Ji^i^i.^{ki, k2, k^)pf\ki, k2, k^)Ai^ {ki)Ai^{k2)Ai^{k3) 



is the CMB bispectrum. The integral, 

can be calculated relatively quickly by using the recurrence relations discussed in the Ap- 
pendix. Eq. (16) provides a general formalism for calculating the bispectrum of the CMB 
starting from any primordial spatial bispectrum. 

FMG defined a statistic Bi which is related to the bispectrum by 

= (2/ +\)3/2 (000) 

and this is especially easy to calculate from Eq. ([16|) . Because of the symmetry, we can 
assume ki < k2 < k^ and introduce three parameters, r, u, and v, defined by, 

kl = r, (19) 
k2 = ur, (20) 
^3 = (m + v)r, (21) 

where u > 1 and < f < 1, to ensure the triangle relation between the three /c's. Finally, 

we get 

('gyj-')3 POO roo 

Bi = Q^^ r^Ai{r)dr Ai{ur) du 

V47r Jo J I 

X / {u + vY dv Jm(^,ur^r{u + v)^Plp'' (^^ur^r{u + v)^Ai{r{u + v)^. (22) 

•J 

The factor 6 comes from permutation of /ci, /c2, and k^. 

At the large angular scales relevant for COBE, the Sachs- Wolfe effect p2[ dominates, 
and we simply have Ai{k) = As^ji{kAr]), where At] = (r^o — 77*) is the conformal time 
between now and the surface of last scatter, and A^^ = 1/3 for a critical-density model with 
primordial adiabatic perturbations. 
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IV. SINGLE-FIELD SLOW-ROLL INFLATION 



(3) 

The formalism outlined above is general and applicable to any form of . However, as 

pointed out by Luo ||2^, the Pjj^^ generated by slow- roll inflation models has a special 
form that simplifies the calculation. In fact, the power-spectrum normalized bispectrum. 



If 



B, 



(23) 



that was extracted from the COBE-DMR data by FMG can be expressed in a very simple 
form for models with the following form of the spatial three-point correlation function: 

Pf k2, h) = fih, k2) + fih, h) + fih, h), (24) 

where f{x, y) is an arbitrary function of x and y. We start with the Sachs- Wolfe formula. 



AT 



fi) = / d^k^{k)e 



iA»7kn 



(25) 



In this case, the angular three-point correlation function becomes 

^(fii, na, lis) = 4w / d^ki d^k2 d^k^ e^^^C^i ^i+k^ ft^+ka-Aa) 

x6d{\^i + k2 + k3)[/(A;i, k2) + fik^, k^) + f{ki, k^)]. 



(26) 



Now, let us calculate the first term associated with f{ki,k2). By integrating out k^, then 
using spherical-harmonic orthonormality, Eqs. (|T3|) and (|1^), and 



^yi™(ni)y^;„(n2) = 5^(111 - fia) 

Im 

we have the bispectrum from the first term of Eq. (|26|) , 

P;,i,i3(firstterm) = AL(47r)4 



(27) 



(2/i + l)(2/2 + l)(2/3 + l) / h k h 

i 



47r 

X j dkidk2klklf{ki,k2)jlikiA7])jl{k2A7]), 

and similarly for the second and third terms. The quantity If then takes the form. 



(28) 



If 



3 / dk, dk2 ki ki f{k,, k2) jf{k, At^) 3t{k2Ar,) 



Jdkk^P;'\k)jf{kAr,) 



3/2 



(29) 



For slow-roll-inflation models, /(fci, ^2) takes the form |[TT| , p!2 



fik^,k2) = A■,^aPi'\k^)Pi'\k2), 



(30) 



where 
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infl 



25mpi 

48 TT 



(31) 



is a constant depending on the height V of the inflaton potential and its first and second 
derivatives, V and V" , respectively, with respect to the infiaton (p |jl2|. Thus, 



If 



3V47rA 



infl 



(32) 



For a flat scale-invariant spectrum of primordial density perturbations, Pjp'\k) = Ank ^, 
where (using our Fourier and scale-factor conventions). 



A 



H 



2h{V'fm\ 



(33) 



PI 



We thus flnd that in terms of the usual slow-roll parameters 



PI 



77 ' ^ 



167r \ V 



our central result for If can be expressed as 





V" 




m 


Svr 


1/ 


-5 





(34) 



(3e-2r7). 



(35) 



PI 



Given that iV^^"^ /to^i) / y/e ~ ~ 10 ^ is the density-perturbation amplitude, and 



e,!] ^ 1 in slow-roll inflation, it is clear that a value + l)/f ~ 1 is inconsistent with 
single-fleld slow-roll inflation. 



V. SINGLE-FIELD MODEL WITH A FEATURE 

Eqs. (|3^) and (|35|) suggest that the bispectrum could be larger if the inflaton potential 
was not as smooth as required for slow-roll inflation. Thus, consider now an inflation model 
with a "feature" in the inflaton potential, such as in those models introduced to account for 
the claimed detection of a feature in the measured mass power spectrum |T^,|TB|,|2^ . More 



precisely, consider a single-fleld inflation model in which the inflaton potential is continuous, 
but its scaled slope 5 = mpi{V'/V) = 4(7re)^/2 jg discont inuous at ko. The derivative of 
the scaled slope (with respect to the fleld 0) is proportional to a Dirac delta function, 
6' = Afcs,tkQTnpl6D{k — ko). The spike occurs at the value that has when the comoving 
wavenumber k^ exits the horizon, and Afeat is the dimensionless amplitude of the delta 
function. 

We must have 5 ^ 1 on both sides of the delta function, and this places a constraint to 
the amplitude Afeat- The change to 5 as the fleld passes through the Dirac delta function is 

= 1 6'd(l) = J ^^^S{k - ko) d<j), (36) 
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and this must be small compared with unity. Using dcf) = {d(j)/ dink) dink ~ {(j)/H){dk/k) 
and ~ {mpy/8TT){V' /V) during slow-roll, we get A5 ~ Afeat^/STr, and so Afc^t ~ Svr. 

To proceed more precisely, we must realize that the derivation |lTT| , p!2| that leads to 
our Eq. (^) is valid only if the inflaton potential varies smoothly. By including Dirac 
delta function in Eq. (32) in Ref. [|1^, we arrive at the correct expression for the spatial 
gravitational-potential bispectrum for an inflaton potential with a feature: 

P!p\ku k2, k,) = ^A_Afc,,A;o5(A;i - k,)P!^\k2)PS\h) + {h ^ k^} + {k, ^ k^}. (37) 

We can still approximate P^"* = Ajjk'^ (since we are still considering only e, -C 1 other 
than the point k = ko), and then from Eqs. ([T6| ) and (plSf), we have 



B, = ^Afcat^i I kldk,kldk2kldk,Jniik,,k,,k3)Mk,)Ai{k2)Mk'^'^^^' ^° 



2,3 l3 
tv2 



^ k2} + {/ci ^ /C3}. (38) 
By the symmetry of ^1,^2, and k^, we obtain 



/ Ai{k2)k^^dk2 I _ Ai{k3)jui{ko,k2,k:i)k^'dks, (39) 



where kx = k^Ar]. Again, during slow roll, the first term H/{m-p\5) ~ H'^/'P ~ 10~^. 
Numerical integration shows that the rest expression peaks at a value of / corresponding 
to the scale of /cq, but with an amplitude ^ O.lAfeat- Therefore, even with a discontinuity 
in the slope of the inflaton potential, the predicted If is still several orders of magnitude 
smaller than unity (the desired value to explain the COBE anomaly found by FMG around 
/ = 16) as long as the slow-roll conditions are satisfied everywhere else. 



VI. DISCUSSION 

We gave a numerically manageable general formalism that allows us to calculate the 
bispectrum of the CMB starting from an arbitrary spatial bispectrum. For certain forms of 
the spatial bispectrum, including those given by single-field slow-roll inflation models, the 
Sachs- Wolfe calculation results in an analytic expression for the non-Gaussian statistic If] 
it turns out to be roughly the density-perturbation amplitude times a linear combination 
of the slow-roll parameters e and r]. This result demonstrates that FMG's non-Gaussian 
signal is inconsistent with single-field slow-roll inflation models. The predicted value of If 
can be increased if the slope of the inflaton potential is discontinuous, but even this larger 
non- Gaussian signal is too small to account for the detection found in Ref. Although 
a complete calculation would require detailed specification of a model, our results for the 
discontinuous single-field model suggest that generic multiple-field models designed to pro- 
duce a break in the slope of the mass power spectrum |lT7| , p!8| will be unable to produce a 
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non-Gaussian signal large enough to explain the COBE result. Counter-examples can be 
found. In particular, inflation models can be designed in which density perturbations are 
produced by quantum fluctuations of some scalar field other than that driving inflation. 



For example, in Ref. ||19[, non-Gaussian isocurvature perturbations to an axion density are 
produced by quantum fluctuations in the axion field, but inflation is driven by some other 
scalar field not associated with Peccei-Quinn symmetry breaking. The mechanism of Ref. 



PT| is somewhat similar. In the model of Ref. |^ , non-Gaussian density perturbations are 
produced (at least in part) by the square of a second scalar field (although details of the 
inflationary dynamics are not presented). In conclusion, if the detection of nonzero If is 
ultimately attributed to the CMB, it raises a serious problem for inflationary models where 
quantum fluctuations in the inflaton give rise to large-scale structure. 
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APPENDIX 

We need to evaluate the following integral: 

Ji^i2i^{ki,k2,k3) = Jhikix) ji2{k2x) ji^{k3x) x"^ dx, (40) 







where ji{x) is spherical Bessel function. The wavenumbers ki, k2, and k^ satisfy the triangle 
relation (they should be able to form a triangle), and li + I2 + h has to be an even number. 
To evaluate the integral, we use the recursive technique discussed in Ref. P3. To begin, we 
evaluate the integral for some special cases: 



00 



T r' ^ 1-,^ \ iuk^x 7 sin A;2a; sin /esx n, 



I even 



-r— — / Piiu)du / e™''''' sin ^23; sin fcsx 

[ Pi{u) du \5{uki + k2 + fcs) - 5{uki - k2 + k^) - 
Jo 



-71 



Ai^k2k3 Jo 
6{uki + k2- k^) + 6{uki - k2 - k^)] 

A relation ^ 



/ Pi{u)du5{uki — \k2 — k^ 
Jo 



Ai^k2k3 Jo 

" ' -P.C^]; (41) 



Ai'- kik2k3 \ ki 
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Jo li'- J-1 k2X k^x 

Zqdd 1 1 /"^ PfoA^o, f°° Jukix 



k2k3 J-1 



/i /-OO 
Pi{u)du I e'''^^'' COS k2X sin k^x dx 
-1 J-00 

[ Pi{u)du\5{uki + k2 + h) + 5{uki - k2 + fcs) - 
Jo 



4ii+^k2k3 

5{uki + k2- ks) - 5{uki - k2 - ks)] 

Are|ation ^ f Pi(u)du[6{uki - /Cs + h) - 6{uki + k2 - h)] 

Ai'-+'-k2ks Jo 
71 1 p f k2 — k^X 



(42) 



To evaluate the integral for other cases, we use the recursion relation, 

21 + 1 

ji-i (a) + ji+i (a) = ji (a) , (44) 

a 

to obtain 



Ai/2/3(^i' ^2, ^3) = i jiAkix) ]i^{k2x) ii.,{k2,x) X dx 

ki 



2/1 + 1 
k2 



(</h-l,/2,/3 + Jll + l,l2,h) (45) 

{Jh,h-i,h + •^h,«2+i,i3)- (46) 



2/2 + 1 

Therefore, we have the recursion relation we need: 

Jh,h+ih = fc^ 2/i + 1 ^"^'^'^''^''3 + Jh+i,i2,h) - Jh,h-i,h- (47) 

If we begin with Eqs. (|41|), (|42D, and (|43|), then we can use Eq. (^T]), to recursively evaluate 
the integral (|40|) for any combination of /i, /2, and /s. 
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